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Abstract.  The  interpretation  of  recent  odd  hydrogen  measurements  in  the  stratosphere 
from  balloons  and  in  the  mesosphere  from  space  indicates  a  serious  lack  of  understanding 
in  atmospheric  HOx  chemistry.  In  order  to  resolve  these  persisting  problems,  coincident 
measurements  of  FIOx  molecules  and/or  measurements  that  cover  both  altitude  regions 
are  desirable.  In  this  work,  the  airborne  2.5  THz  heterodyne  spectrometer  Terahertz  OH 
Measurement  Airborne  Sounder  (THOMAS)  is  introduced.  Since  the  first  THOMAS 
measurements  in  1994/1995,  the  spectrometer  was  significantly  improved  by  modification 
or  replacement  of  individual  components.  The  THOMAS  instrumental  setup  and  properties 
are  presented  together  with  a  retrieval  algorithm  for  atmospheric  parameters  based  on  a 
Phillips-Tikhonov  regularization  scheme.  Furthermore,  the  results  of  a  complete  error 
assessment  are  given.  In  August  1997,  during  the  second  CRISTA/M AHRSI  campaign 
(Cryogenic  Infrared  Spectrometers  and  Telescopes  for  the  Atmosphere  /  Middle  Atmo¬ 
sphere  High  Resolution  Spectrograph  Investigation),  OH  observations  were  performed 
by  THOMAS  covering  altitudes  between  about  30  and  90  km  over  a  full  diurnal  cycle. 
Hydroxyl  column  densities  derived  from  THOMAS  measurements  are  presented  and  com¬ 
pared  to  photochemical  model  results.  The  model  calculations  using  the  standard  HOx 
chemistry  systematically  show  higher  values  by  about  1 5  %  for  the  40-90  km  and  50-90  km 
OH  columns.  Moreover,  a  recently  proposed  change  of  an  HOx  chemistry  reaction  rate 
is  included  into  the  comparison  which,  for  the  same  altitude  intervals,  yields  OH  column 
densities  that  are  about  10%  lower  than  the  THOMAS  measurements.  A  detailed  compar¬ 
ison  of  the  THOMAS  and  MAHRSI  measurements  is  presented  in  a  seperate  publication 
[Englert  et  al.,  2000], 


1.  Introduction 

Ozone  chemistry  throughout  the  middle  atmosphere 
(stratosphere  and  mesosphere)  is  closely  connected  to  the 
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odd  hydrogen  family  (HOx  -  H,  OH,  HO2).  A  detailed  un¬ 
derstanding  of  atmospheric  chemistry  involving  odd  hydro¬ 
gen  is  therefore  of  great  importance  for  the  refinement  of 
atmospheric  models  that  enable  us  to  understand  the  current 
state  of  the  atmosphere  and  possible  future  changes.  Com 
paring  photochemical  model  results  to  measurements  is  a 
rigorous  test  of  the  current  understanding  of  atmospheric 
chemistry.  However,  measuring  these  radicals  with  suffi¬ 
cient  precision  and  accuracy  poses  an  experimental  chal¬ 
lenge.  Owing  to  major  advances  in  the  development  of  suit¬ 
able  in  situ  and  remote  sensing  techniques,  a  number  of  in¬ 
struments  have  supplied  valuable  data  in  the  last  decades. 

Balloon-borne  measurements  of  OH  concentration  pro¬ 
files  in  the  middle  and  upper  stratosphere  have  been  per¬ 
formed  using  Lidar  [Heaps  and  McGee.  1983.  1985],  in¬ 
frared  Fourier  transform  spectrometry  f  Carli  and  Park,  1988; 
Carli  et  al..  1989;  Park  and  Carli,  1991;  Johnson  et  al., 
1995],  ultraviolet  spectrometry  [Torr  et  al.,  1987],  and  in¬ 
frared  Fabry-Perot  spectrometry  [Pickett  and  Peterson, 
1993],  A  comparison  of  some  results  can  be  found  [e.g., 
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Pickett  and  Peterson,  1993],  Time-dependant  HO;c  profiles 
of  the  middle  stratosphere  resulting  from  a  balloon-borne  in¬ 
frared  spectrometer  experiment  were  presented  for  the  first 
time  by  Chance  et  al.  [1996].  The  HO^  distribution  was 
found  to  largely  agree  with  the  standard  photochemical  mod¬ 
el  except  for  OH  at  the  highest  measurement  altitudes  near 
38  km  which  was  significantly  higher  («  33%  or  3<r  mea¬ 
surement  uncertainty)  than  the  model  values.  (Descriptions 
of  photochemical  models  are  given  by  Siskind  et  al.  [1995], 
Pickett  and  Peterson  [1996],  Salawitch  et  al.  [1994],  Clan¬ 
cy  et  al.  [1994J,  or  Sandor  and  Clancy  [1998].)  Pickett 
and  Peterson  [1996]  reported  stratospheric  OH  concentra¬ 
tion  profiles  retrieved  from  balloon-borne  far-infrared  limb 
observations  of  a  triple  Fabry-Pcrot  interferometer.  The  ob¬ 
servations  cover  full  diurnal  cycles  and  basically  agree  well 
with  the  applied  photochemical  model.  However,  the  com¬ 
parison  of  a  data  set  from  Pickett  and  Peterson  [  1 996]  with 
a  different  model  shows  an  OH  underestimation  of  about 
20  %  (w  2  cr  measurement  uncertainty)  at  37  km  [Ostermann 
et  al.,  1997],  Long-term  ground-based  total  column  mea¬ 
surements  have  been  reported  by  Burnett  and  Minschwaner 
[1998J.  They  observe  significant  changes  in  total  column 
OH  during  the  last  2  decades.  Since  tropospheric  OH  is  only 
a  minor  contribution  to  the  total  column,  they  conclude  that 
these  changes  are  most  likely  originating  from  changes  in 
stratospheric  or  mesospheric  OH.  However,  the  measured 
OH  column  abundance  is  significantly  higher  than  that  pre¬ 
dicted  from  standard  chemistry. 

Hydroxyl  measurements  in  the  mesosphere  from  a  space 
platform  have  been  performed  with  the  Middle  Atmosphere 
High  Resolution  Spectrograph  Investigation  (MAHRSI)  in 
1994  [Conway  et  al,  1996]  and  1997.  The  reported  OH  con¬ 
centration  profiles  are  inconsistent  with  model  calculations 
using  the  standard  photochemistry  (later  referred  to  as  mod¬ 
el  A).  Two  alternative  changes  to  standard  chemistry  were 
proposed  to  explain  the  MAHRSI  observations  [Summers  et 
al.,  1997],  Either  a  50%  reduction  of  the  ki  rate  coefficient 
for 

0  +  H02^0H  +  02  (1) 

(later  referred  to  as  model  B)  or  a  20  %  reduction  in  the  rate 
coefficient  ki  together  with  a  30  %  increase  in  the  fc2  rate 
coefficient  for 

0H  +  H02-tH20+O2  (2) 

lead  to  consistency  of  OH  model  results  and  observations. 
Compared  to  the  standard  model,  this  rate  constant  reduc¬ 
tion  also  resulted  in  a  substantially  better  agreement  between 
the  model  ozone  concentration  and  simultaneous  ozone  mea¬ 
surements  at  50  km  by  Cryogenic  Infrared  Spectrometers 
and  Telescopes  for  the  Atmosphere  (CRISTA).  Ground- 
based  observations  of  lower  mesospheric  H02  [Sandor  and 
Clancy,  1998]  favor  the  reduction  of  k\  only.  This  modi¬ 
fication  of  the  HOj,  chemistry  has  only  a  minor  effect  on 
lower  stratospheric  OH,  H02,  and  O3  abundances  [Sandor 
and  Clancy,  1998;  Summers  et  al,,  1997;  Chen  et  al.,  1997], 
However,  a  change  of  k\  by  50%  would  mean  a  change  of 
more  than  two  0  measurement  uncertainty  and  has  been  op¬ 


posed  by  the  laboratory  and  stratospheric  modeling  commu¬ 
nity  [cf.  Sandor  and  Clancy,  1998].  The  MAHRSI  data 
also  corroborate  the  existance  of  an  unexplained  H20  layer 
between  65  km  and  70  km  observed  by  the  Halogen  Occul- 
tation  Experiment  (HALOE)  [Summers  et  al.,  1997;  Siskind 
and  Summers,  1998]. 

Simultaneous  measurements  of  OH.  H02,  H20,  and  O3  in 
the  upper  stratosphere  have  recently  been  reported  by  Jacks 
et  al.  [1998].  Above  35  km  they  found  OH  concentrations 
in  agreement  with  standard  chemistry,  while  the  observed 
H02  concentrations  showed  about  25  %  higher  values  than 
the  standard  model.  The  enhanced  H02  abundance  seems 
to  be  consistent  with  the  mesospheric  H02  observations  of 
Sandor  and  Clancy  [1998].  The  measurements  reported  by 
Jucks  et  al.  [19981  imply  a  higher  total  HO*  concentration 
as  well  as  a  difference  in  the  partitioning  of  OH  and  H02. 
In  conclusion,  the  present  set  of  HO.r  observations  neither 
shows  general  mutual  agreement  nor  good  agreement  with 
standard  model  calculations. 

In  August  1997,  during  the  second  MAHRSI  mission,  the 
2.5  THz  OH  Measurement  Airborne  Sounder  (THOMAS) 
performed  vertical  sounding  measurements  of  OH  between 
about  30  km  and  90  km.  Tire  geometry  of  the  observations 
was  chosen  in  a  way  that  they  coincided  in  space  and  time 
with  MAHRSI  measurements  whenever  possible,  thus  al¬ 
lowing  a  direct  intercomparison.  Together  with  THOMAS, 
the  Fast  In  Situ  Stratospheric  Hygrometer  (FISH)  [Zdger, 
1996]  was  operated  on  board  the  Deutschen  Zentrum  fiir 
Lult-  und  Raumfahrt  (DLR)  research  aircraft  FALCON  and 
supplied  simultaneous  water  vapor  data  at  flight  altitude  for 
most  of  the  THOMAS  observations. 

The  next  section  of  this  paper  presents  some  fundamen¬ 
tals  of  far-infrared  spectroscopy.  The  THOMAS  instrument 
and  the  measurement  conditions  are  described  in  the  follow¬ 
ing  section  with  emphasis  on  the  major  improvements  of  the 
instrument  made  after  the  first  flights  of  THOMAS  in  1994 
[77ft  et  al.,  I995a[  and  1995  [77ft  et  al.,  1995b],  The  in¬ 
version  algorithm  is  outlined  in  section  4  together  with  an 
assessment  of  possible  errors  in  the  experimental  data  as 
well  as  in  the  inversion  of  OH  concentration  profiles  from 
the  measured  far-infrared  emission  spectra.  In  sections  5 
and  6  the  inversion  results  are  shown  and  compared  to  mod¬ 
el  calculations,  A  detailed  comparison  to  results  from  the 
second  MAHRSI  campaign  is  presented  in  a  seperate  publi¬ 
cation  [Englert  et  al,  2000]. 

2.  Far-Infrared  Thermal  Emission 
of  Atmospheric  OH 

In  a  feasibility  study  by  Miller  et  al.  [1992]  the  pure  rota¬ 
tional  transitions  of  OH  at  2.5  THz  (83  cm-1)  were  found 
to  be  suitable  for  measuring  vertical  OH  profiles  with  an 
airborne,  uplooking  heterodyne  receiver  flying  at  altitudes 
above  the  tropopausc  in  order  to  reduce  signal  attenuation 
by  water  vapor  absorption.  The  double  side  band  receiver 
THOMAS  observes  these  thermal  emission  lines  of  OH  in 
the  lower  side  band  superimposed  with  an  H20  transition  at 
84.456  cm-1  in  the  upper  side  band  that  provides  informa¬ 
tion  about  the  H20  abundance. 
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The  radiance  I  at  wavenumber  v  received  by  THOMAS 
is  described  by  the  radiative  transfer  equation.  Its  integral 
form,  neglecting  scattering  and  assuming  local  thermody¬ 
namic  equilibrium,  is  given  by  [Liou,  19801 


=  /<*»  7>,oo) 


-  J~  dz' B[v,T(z')]  dT{Qzf]  ,  (3) 


where  B  is  the  Planck  function  at  temperature  T  and  is  a 
background  contribution,  for  example,  the  3  K  background 
radiation  of  deep  space.  The  monochromatic  transmission 
T  (relative  to  the  observer  at  z )  is  given  according  to  Beer’s 
law  by 


T(y,z’)  =  expj ~  (E  krm(v,z")nm(z") 

+  acm{i/,  z")^  dz"  j,(4) 

where  nm  is  the  number  density  of  molecule  m  and  is  its 
absorption  cross  section  for  resonant  absorption.  The  con¬ 
tinuum  absorption  due  to  nonresonant  effects  is  denoted  by 
the  continuum  absorption  coefficient  acm.  Note  that  the  ab¬ 
sorption  cross  section  and  the  absorption  coefficient  depend 
on  pressure  and  temperature,  but  for  brevity  we  use  the  con¬ 
densed  notation  fcjJ,(i/,2)  =  k^\y,p{z),T(z)\.  For  sim¬ 
plicity,  we  have  assumed  a  vertical  path  geometry,  but  (3) 
and  (4)  are  easily  generalized  to  arbitrary  slant  paths.  Instru¬ 
ment  influences  on  the  measured  spectrum  are  not  consid¬ 
ered  in  equation  (3). 


2.1.  Resonant  and  Nonresonant  Absorption 


For  an  individual  line  l  the  spectral  absorption  cross  sec¬ 
tion  is  the  product  of  the  temperature-dependent  line 
strength  Smj(T)  and  a  normalized  line  shape  function  b  de¬ 
scribing  the  broadening  mechanism,  krml{v,z)  = 


Smj(T)b  v,i'm,h7m,l(p,  r)J,  where  i/m,j  is  the  line  cen¬ 


ter  and  7m,i  is  the  line  width.  In  general,  the  absorption 
cross  section  is  obtained  by  summing  over  the  contributions 
from  many  lines,  fcJJ,  =  52j  krm  l.  For  the  far-infrared  and 
under  atmospheric  conditions,  the  combined  effect  of  pres¬ 
sure  broadening  (corresponding  to  a  Lorentzian  line  shape) 
and  Doppler  broadening  (corresponding  to  a  Gaussian  line 
shape)  can  be  represented  by  a  Voigt  line  profile. 

Nonresonant  absorption  is  described  by  an  absorption  co¬ 
efficient  acm  showing  a  weak  wavenumber  dependence.  In 
the  spectral  region  observed  by  THOMAS  the  H2  0  continu¬ 
um  and  the  N2  continuum  are  dominant  [ Biihler  et  al. ,  1 996; 
Urban,  1998],  The  semiempirical  H20  continuum  described 
by  Clough  et  al.  [1989]  is  a  pure  far  wing  continuum  cor¬ 
recting  the  contribution  of  lines  more  than  25  cm-1  away 
from  the  spectral  region  of  interest.  Liebe  et  al.  [  1 993]  for¬ 
mulated  an  empirical  expression  for  the  continuum  contri¬ 
bution  of  collision  induced  N2  absorption.  Because  of  the 
lack  of  experimental  data  in  the  spectral  region  of  interest, 
modeling  these  continua  is  connected  with  high  uncertain¬ 
ties.  Especially  the  N2  continuum  is  less  established  than  the 


Figure  1.  Comparison  of  H20  and  N2  continuum  absorption 
coefficients 

H20  continuum,  but  its  application  yields  better  agreement 
between  measured  spectra  and  radiative  tranfer  calculations. 
Figure  1  shows  both  continua  calculated  for  the  Air  Force 
Geophysics  Laboratory  (AFGL)  U.S.  standard  atmosphere 
[Anderson  et  al.,  1986].  In  wide  regions  of  the  atmosphere 
the  N2  continuum  is  dominant  except  for  altitudes  near  the 
flight  altitude  of  THOMAS,  where  both  continua  are  sim¬ 
ilar  in  magnitude.  The  H>0  continuum  is  proportional  to 
the  II20  concentration.  At  flight  altitudes  of  the  THOMAS 
instrument  near  the  tropopause,  where  the  H20  concentra¬ 
tion  is  exponentially  decreasing  with  increasing  altitude,  the 
H20  concentration  is  highly  variable,  and  the  H20  contin¬ 
uum  becomes  dominant  in  regions  of  high  H20  concentra¬ 
tion.  The  main  influence  of  the  continua  on  the  OH  lines  is 
due  to  damping  of  the  signal  by  the  continuum  absorption 
on  the  first  few  kilometers  above  the  observer  height. 

2.2.  Spectral  Information  Contents  for  Retrieval 

Spectral  information  about  the  OH  concentration  profile 
can  be  found  in  the  line  shape  of  the  OH  lines.  For  altitudes 
lower  than  50  km,  where  pressure  broadening  is  dominant 
for  OH,  the  line  width  of  the  emission  lines  decreases  with 
increasing  altitude.  The  pressure  broadening  line  width  is 
proportional  to  the  pressure  and  ranges  from  2x  10-2  cm-1 
at  flight  altitude  to  lxl0~4  cm-1  at  50  km.  As  a  conse¬ 
quence,  the  line  profile  has  to  be  completely  resolved  to  al¬ 
low  the  retrieval  of  altitude  information.  Above  50  km  the 
nearly  altitude-independent  Doppler  broadening  becomes 
dominant  so  that  the  line  shape  of  the  OH  lines  contains  in¬ 
formation  about  the  mesospheric  column  only. 

3.  Improved  THOMAS  Instrument 

The  THOMAS  instrument  is  a  further  development  of  a 
2.5  THz  heterodyne  spectrometer  of  the  Max-Planck-Insti- 
tute  for  Radioastronomy  in  Bonn,  Germany,  where  it  was 
built  for  astronomical  research  on  board  the  Kuiper  Air¬ 
borne  Observatory  [ Roser ,  1991].  THOMAS  is  a  similar 
spectrometer,  modified  for  atmospheric  measurements  using 
the  DLR  research  aircraft  FALCON.  After  two  campaigns  in 
1994  and  1995,  major  components  of  THOMAS  have  been 
optimized  or  replaced  to  improve  the  spectrometers  perfor¬ 
mance. 
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Figure  2  shows  a  schematic  of  the  THOMAS  instrument. 
The  atmospheric  signal  enters  the  fuselage  through  a  spe¬ 
cial  aircraft  window.  After  being  reflected  by  the  pointing 
mirror  and  a  plane  mirror,  one  half  of  the  unpolarized  radi¬ 
ation  gets  reflected  into  the  diplexer  by  the  90°  wire  grid. 
About  90  %  of  the  monochromatic  radiation  of  the  first  local 
oscillator  (LO),  a  CO2  laser  pumped  methanol  laser  (i'lq 
=  84.15095cm-1)  is  transmitted  by  the  90°  wire  grid  into 
the  diplexer.  The  remaining  fraction  is  reflected  and  detect¬ 
ed  by  a  chopped  pyroelectrical  detector  to  monitor  the  LO 
performance.  Tire  diplexer  rotates  the  polarization  of  the 
first  LO  to  the  polarization  of  the  atmospheric  signal.  The 
combination  of  the  first  LO  beam  and  the  atmospheric  radi¬ 
ation  is  focused  onto  an  open  structure  corner  cube  mixer  by 


an  off-axis  parabolic  mirror.  The  actual  mixer  consists  of 
a  whisker  antenna  contacting  a  honeycomb  Schottky  diode. 
The  diode  type  T23,  supplied  by  the  University  of  Virginia, 
was  used  for  this  work  [Crowe  et  al.,  1992],  The  resulting 
double  side  band  (DSB)  spectrum  around  8.5  GHz  is  am¬ 
plified  and  mixed  down  to  a  center  frequency  of  2.1  GHz 
(bandwidth:  1.4  GHz)  by  the  intermediate  frequency  chain 
(IF  chain).  An  acousto  optical  spectrometer  (AOS)  is  used  to 
perform  a  spectral  analysis  of  the  signal.  The  resulting  spec¬ 
tra  are  displayed  and  stored  by  a  personal  computer.  The 
computer  also  stores  data  from  the  aircraft  interface  like  al¬ 
titude,  roll  and  pitch  angle,  outside  temperature,  and  head¬ 
ing  as  well  as  instrumental  housekeeping  data  like  CO2  laser 
power,  methanol  laser  power,  mixer  voltage,  temperatures  of 
the  window,  the  hot  and  cold  reference  black  bodies,  and  the 
total  power  leaving  the  IF  chain.  The  radiometric  calibration 
of  the  spectra  is  performed  using  spectra  of  the  two  reference 
black  bodies  at  liquid  nitrogen  and  ambient  temperature. 

Critical  components  and  parameters  of  THOMAS  instru¬ 
ment  are  discussed  in  the  following,  including  the  modifica¬ 
tions  that  were  made  since  the  1994/1995  flights. 

3.1.  Quasi-Optical  Design 

The  quasi  optics  of  the  front  end  of  THOMAS  has  been 
designed  for  a  beamwaist  of  9  mm  corresponding  to  a  Ray¬ 
leigh  length  of  2.157  m  and  an  asymptotic  field  of  view  angle 
of  0.24°.  The  defacto  field  of  view  depends  on  the  whisker 
position  in  the  corner  cube  structure  of  the  mixer,  the  po¬ 
sition  of  the  mixer  relative  to  the  parabolic  mirror,  and  the 
qualitiy  of  the  wire  grids  (compare  with  section  3.3).  For  the 
atmospheric  measurements  the  mixer  position  was  adjusted 
for  maximum  reception  of  the  first  LO.  An  adjustment  of 
the  whisker  position  was  not  possible  during  operation  of 
the  mixer.  Since  the  field  of  view  is  not  as  critical  for  a  ver¬ 
tical  sounding  instrument  like  THOMAS  compared  to  limb 
sounding  experiments,  this  adjustment  procedure  was  con¬ 
sidered  to  be  sufficient  (see  also  section  4.2).  Because  of  the 
narrow  field  of  view  in  combination  with  the  uplooking  ge¬ 
ometry,  a  pencil  beam  was  assumed  for  the  retrieval  of  OH 
concentrations. 

3.2.  Aircraft  Window 

For  the  earlier  atmospheric  measurements  with 
THOMAS  a  polyethylene  (PE)  aircraft  window  had  been 
used  [Titz  etal. ,  1995a],  The  transmittance  of  the  PE  aircraft 
window  was  only  61  %,  and  a  proper  temperature  measure¬ 
ment  of  the  window  was  difficult  because  of  the  low  ther¬ 
mal  conductivity  of  PE  that  caused  large  temperature  gradi¬ 
ents  in  the  window  during  the  flight.  Because  of  these  dis¬ 
advantages,  the  polyethylene  aircraft  window  was  replaced 
by  an  antireflection-coated,  wedged,  single-crystal  silicon 
window  [Englert  et  al. ,  1999a].  The  free  aperture  of  the 
new  window  is  5  cm  x  15  cm  allowing  pointing  angles  be¬ 
tween  5°  and  20°  toward  horizontal.  The  temperature  of 
the  new  window,  which  is  needed  to  account  for  the  far- 
infrared  radiation  emitted  by  the  window,  can  be  measured 
more  accurately  due  to  the  much  higher  thermal  conductivi- 
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ty  of  silicon.  Moreover,  the  absorption  coefficient  of  silicon 
in  the  far-infrared  is  significantly  lower,  making  the  window 
more  transparent.  Laboratory  measurements  showed  that  the 
transmittance  of  new  window  is  about  89  %.  The  losses  can 
be  divided  into  about  7  %  reflection  losses  and  about  4  % 
losses  due  to  absorption.  However,  the  transmittance  of  the 
window  decreased  up  to  about  20  %  while  it  was  illuminated 
by  direct  sunlight  during  the  measurement  flights.  Temper¬ 
ature  changes  of  the  window,  typically  of  the  order  of  2  K, 
can  be  ruled  out  as  a  reason  for  the  transmittance  change. 
Laboratory  measurements  with  THOMAS  and  with  a  Fouri¬ 
er  transform  spectrometer  [Englert  et  al.,  1999b]  confirmed 
that  the  change  in  absorption  is  due  to  an  increased  number 
of  electrons  in  the  semiconductor’s  conduction  band.  The 
source  of  the  additional  free  charges  in  the  silicon  are  tran¬ 
sitions  of  electrons  from  the  valence  band  to  the  conduc¬ 
tion  band  caused  by  radiation  with  energies  greater  than  the 
bandgap  (visible  light,  UV). 

A  correction  of  the  photon  induced  absorption  was  de¬ 
termined  using  two  successive  airborne  measurements,  one 
with  and  one  without  direct  sunlight  illuminating  the  win¬ 
dow.  These  particular  measurements  have  been  selected 
from  all  available  spectra  because  the  coinciding  in  situ  wa¬ 
ter  vapor  measurements  of  the  FISH  instrument  aboard  the 
FALCON  showed  that  the  flight  altitude  was  well  above  the 
tropopause  for  this  day.  Therefore  the  effect  on  the  measured 
spectrum  due  to  changing  water  vapor  abundance  in  the  first 
few  kilometers  of  the  line  of  sight  is  expected  to  be  mini¬ 
mal,  and  the  observed  change  in  signal  absorption  can  be  at¬ 
tributed  to  photon  induced  absorption  only.  The  absorption 
change  of  the  window  was  basically  determined  by  fitting 
the  window  absorption  so  that  the  shapes  of  the  measured 
water  vapor  emission  lines  matched  best.  This  procedure 
assumes  identical  distribution  of  middle  atmospheric  water 
vapor  for  both  measurements. 

Figure  3  shows  the  two  spectra  measured  with  and  without 
the  Sun  shining  directly  onto  the  window.  In  the  upper  panel 
of  Figure  3  the  raw  spectra  are  displayed,  whereas  in  the 
lower  panel  the  window  corrected  spectra  are  shown.  The 
predominant  features  in  the  spectra  are  a  rotational  line  of 
water  vapor  at  about  84.456  cm-1  in  the  upper  side  band 
and  a  rotational  line  triplet  of  hydroxyl  around  83. 869  cm-1 
in  the  lower  side  band. 

All  the  measured  spectra  have  been  window  corrected  un¬ 
der  the  assumption  of  a  linear  relation  between  the  additional 
absorption  and  the  solar  energy  flux  in  the  visible  and  UV  at 
the  window.  The  diffuse  solar  energy  flux  in  the  shade  was 
estimated  to  be  1 0  %  of  the  direct  radiation  of  the  Sun. 

3.3.  Wire  Grids 

Laboratory  measurements  of  the  THOMAS  far-field  an¬ 
tenna  pattern  showed  that  it  suffers  significant  distortion 
when  reflected  by  the  90°  wire  grid  (wire  diameter:  10  /rm; 
specified  wire  spacing:  d= 25  gm).  This  effect  was  due  to 
large  tolerances  in  the  spacing  of  the  grid  (0  fim  <  d  < 
80  /um).  The  tungsten  wire  grids  have  been  replaced  by  new 
grids  with  a  smaller  tolerance  in  the  grid  spacing  (0  pro  < 
d  <  40  pm),  thus  improving  the  beam  pattern. 


3.4.  First  Local  Oscillator 

Power  stability  and  especially  frequency  stability  is  need¬ 
ed  for  the  first  LO.  The  stability  of  the  FIR  laser  power  is 
maintained  by  the  electronic  adjustment  of  the  CO2  laser  res¬ 
onator  length,  keeping  the  diode  voltage,  which  is  depending 
on  the  LO  power,  at  a  constant  value.  No  frequency  stabi¬ 
lization  of  the  first  LO  is  used.  Measurements  on  board  the 
aircraft  as  well  as  in  the  laboratory  showed  that  the  max¬ 
imum  change  in  LO  frequency,  mostly  due  to  thermal  ex¬ 
pansion  of  the  FIR  laser  resonator,  is  of  the  order  of  2  MHz 
(«  6.7  X 10~°  cm-1).  Note  that  a  drift  in  frequency  of  the 
first  LO  causes  a  broadening  of  the  instrumental  line  shape 
function  (ILS). 

3.5.  Diplexer  Adjustment 

A  position  adjustment  of  one  roof  top  mirror  of  the  diplex¬ 
er  was  performed  immediately  after  takeoff  for  every  flight. 
The  adjustment  was  done  by  maximizing  the  diode  voltage, 
thus  maximizing  the  transmittance  of  the  first  LO.  The  max¬ 
imum  uncertainty  of  the  adjustment  was  of  the  order  of  2  % 
in  diode  voltage  resulting  in  a  sideband  ratio  of  1.00  ±  0.08. 
For  a  misaligned  diplexer  the  sideband  ratio  remains  unity  at 
the  center  frequency  of  the  spectrum  and  increases/decreases 
toward  the  lower/upper  end  of  the  spectrum.  Sideband  ratios 
inferred  from  laboratory  measurements  of  methanol  emis¬ 
sion  lines  gave  no  indication  of  any  major  additional  side¬ 
band  ratio  variations,  for  example,  due  to  changing  mixer 
properties. 

3.6.  First  Mixer  and  Low-Noise  Amplifier 

A  quasi-optical  mixer  structure  introduced  by  Roser 
[1991]  was  used  for  the  measurements  reported  here.  The 
liquid  nitrogen  cooled  low  noise  amplifier  (LNA)  was  re- 
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Figure  3.  (top)  Raw,  radiometric  calibrated  spectra.  The 
dotted  spectrum  was  measured  while  the  Sun  was  shining 
directly  onto  the  aircraft  window.  The  solid  line  spectrum 
was  observed  while  the  aircraft  window  was  in  the  shade, 
(bottom)  Window  corrected  spectra. 
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placed  by  a  significantly  smaller  amplifier  unit  that  can  be 
operated  at  room  temperature  without  any  performance  re¬ 
duction  of  the  spectrometer  [Nitsche,  1998]. 

3.7.  Intermediate  Frequency  Chain 

The  IF  chain  consists  of  standard  microwave  components 
except  for  the  LNA  mentioned  above.  Figure  2  shows  the 
individual  devices.  After  the  second  mixer  stage,  a  power 
meter  is  used  to  measure  the  relative  power  entering  the  AOS 
(middle  frequency:  2.1  GHz,  bandwidth:  1.4  GHz)  during  all 
the  atmospheric  and  reference  black  body  observations.  Tire 
noise  temperature  of  the  IF  chain  was  144  K. 

3.8.  Acousto  Optical  Spectrometer 

Changes  in  the  optics  of  the  AOS  and  a  new  readout  and 
data  processing  unit  significantly  improved  the  performance 
of  the  AOS  that  was  used  together  with  THOMAS  in  the 
1995  campaign.  Subsequent  spectral  averaging  over  three 
CCD  pixels  was  used  to  minimize  the  influence  of  speck¬ 
les  and  to  increase  the  signal  to  noise  ratio.  The  resulting 
full  halfwidth  of  the  instrumental  line  shape  function  was 
reduced  to  about  3  MHz.  The  influence  of  a  possible  fre¬ 
quency  drift  of  the  first  LO  by  about  2  MHz  (compare  with 
section  3.4)  translates  to  a  variation  of  the  half  width  half 
maximum  of  the  final  ILS  by  about  0.2  MHz  (compare  with 
Tablet.).  The  spectral  variation  of  the  ILS  was  neglected 
regarding  the  influence  of  the  first  LO  on  the  ILS. 

3.9.  Calibration 

The  atmospheric  spectral  radiance  cannot  be  measured 
directly  by  submillimcter  spectrometers.  Instead,  a  spec¬ 


tral  noise  equivalent  power  (NEP)  is  measured  which  is  the 
sum  of  the  NEPrcc  (receiver  noise  equivalent  power)  and 
the  NEPsig  (signal  noise  equivalent  power).  The  atmospher¬ 
ic  NEP  measurement  has  to  be  calibrated  to  obtain  the  de¬ 
sired  spectral  radiance.  This  can  be  done  using  measure¬ 
ments  of  two  reference  blackbodies  with  known  tempera¬ 
ture.  The  calibration  procedure  makes  use  of  the  constancy 
of  NEPrec  as  well  as  the  linear  relation  between  NEPsig  and 
the  detected  radiance.  (The  linearity  of  the  instrument  re¬ 
sponse  was  investigated  using  laboratory  measurements  of 
methanol  emission  lines.)  In  other  words,  the  stability  of 
the  instrument  is  needed.  To  obtain  a  minimum  uncertain¬ 
ty  in  the  calibrated  spectrum,  the  temperatures  of  the  ref¬ 
erence  blackbodies  have  to  be  choosen  in  a  way  that  the 
unknown  (atmospheric)  radiance  lies  between  the  radiances 
emitted  by  the  two  blackbodies.  The  optimal  ratio  of  integra¬ 
tion  intervals  for  the  three  individual  measurements  depends 
on  the  blackbody  temperatures  and  the  atmospheric  radiance 
[Klein,  1993].  It  can  be  shown  that  for  the  above  mentioned 
conditions,  the  signal-to-noise  ratio  of  the  calibrated  spec¬ 
trum  is  independent  of  the  blackbody  temperatures.  For  the 
measurements  presented  in  this  work  the  standard  measure¬ 
ment  sequence  (hot-atmosphere-cold-atmosphere. . . )  has 
been  chosen  with  the  following  integration  intervals:  fhot  = 
1 .5  s,  fat m  =  2.0  s,  and  f0Otd  =  2.5  s.  These  time  periods  have 
been  selected  under  consideration  of  the  system  stability  and 
the  time  required  to  turn  the  pointing  mirror.  The  system  sta¬ 
bility  (compare  with  section  3.10)  asks  for  short  integration 
times  so  that  instrument  drifts  can  be  neglected.  On  the  other 
hand,  short  integration  intervals  require  frequent  turning  of 
the  pointing  mirror  which  consumes  valuable  measurement 


Table  1.  Summary  of  Error  Sources  and  Their  Effect  on  the  Weighted  OH  Column  Densities 


Parameter 

Estimated 

Error  in  Weighted  OH  Column 

Magnitude 
of  Variation 

30-50  km 

50-90  km 

40-90  km 

Gaussian  noise  (s.  d.) 

0.75  K 

1.6xl012  cm-2 

1.6xl012  cm-2 

2.3  xl012cm~2 

Regularization  parameter  Aqh 

+100%/- 50% 

15.0% 

5.7% 

6.6% 

ch2o/coh  ratio 

+100%/-50% 

13.0% 

4.4% 

1.4% 

Temperature  profile  bias 

5K 

1.1% 

3.5% 

3.1% 

Ozone  profile 

20% 

2.8% 

<  1% 

<1% 

Pressure  profile  (altitude  shift) 

13% 

8.0% 

4.4% 

4.0% 

H20  variation  at  observer 

10->30ppm 

2.0% 

<1% 

0.3% 

N2  continuum  absorption 

50% 

4.2% 

2.2% 

1.4% 

HjO  continuum  absorption 

50% 

2.8% 

1.3% 

0.9% 

H20  pressure  broadening 
parameter 

15% 

4.3% 

<1% 

<1% 

OH  triplet  line  strength 

1% 

1.0% 

1 .0% 

1.0% 

OH  pressure  broadening 
parameter 

4% 

<1% 

1.3% 

1.1% 

OH  pressure  broadening 
exponent 

20% 

0.4% 

0.9% 

1.1% 

Spectrum  baseline  offset 

3  K 

5.6% 

6.6% 

5.3% 

Sideband  ratio 

8% 

3.0% 

5.6% 

4.6% 

Window  transmittance 

2% 

4.9% 

2.1% 

1.7% 

ILS  width 

0.2  MHz 

1.0% 

0.9% 

0.3% 

AOS  nonlinearity  correction 

0.1 

3.5% 

4.0% 

4.0% 

Pointing  bias 

1° 

5.4% 

7.0% 

6.3% 

Root  mean  square  of  relative 
and  systematic  errors 

~24% 

a  15% 

«14% 
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time.  For  this  work,  the  spectra  have  been  calibrated  with 
a  procedure  similar  to  the  “weighted  calibration”  reported 
by  Crewell  [1993],  This  method  uses  two  hot  and  cold  ob¬ 
servations  (one  preceding  and  one  following  the  atm  mea¬ 
surement)  for  the  calibration  of  a  single  atmospheric  spec¬ 
trum,  thus  canceling  out  linear  instrumental  drifts.  Consecu¬ 
tive  calibrated  spectra  have  been  averaged  to  obtain  a  higher 
signal-to-noise  ratio.  Here  the  individual  spectra  have  been 
averaged  to  give  a  total  measurement  time  of  20  min  or  more 
for  a  final  spectrum. 

3.10.  System  Stability 

Stability  is  a  key  property  of  a  heterodyne  spectrometer 
because  it  is  required  to  determine  the  timing  of  the  mea¬ 
surements.  As  mentioned  above,  three  measurements  are 
needed  to  derive  a  calibrated  spectrum  (compare  with  sec¬ 
tion  3.9).  During  the  time  the  three  data  sets  are  acquired, 
total  system  stability  is  desired.  The  Allan  variance  is  com¬ 
monly  applied  to  investigate  the  stability  of  submillimeter 
spectrometers  [Stanley,  1994;  Allan  1966].  Allan  variances 
have  been  calculated  for  a  THOMAS  airborne  measurement 
using  only  the  observations  of  the  cold  reference  blackbody. 
The  total  measurement  time  was  about  1  hour.  The  results 
for  three  spectral  channels  corresponding  to  pixels  of  the 
AOS  charge-coupled  device  (CCD)  and  for  the  integral  total 
power  measurements  are  depicted  in  Figure  4. 

Minima  of  the  Allan  variance  indicate  the  number  of  mea¬ 
surements  that  can  be  averaged  without  significant  influence 
of  instrument  drifts  on  the  averaged  measured  value.  Said 
differently,  they  determine  the  maximum  interval  of  stabil¬ 
ity.  Figure  4  shows  that  for  this  measurement  stability  was 
accomplished  for  a  longer  time  than  the  period  between  two 


1  10 

Number  of  Averaged  Measurements 


Figure  4.  Allan  variances  for  three  CCD  pixels  (pixel  500, 
close  to  water  vapor  line  center;  pixel  1000,  center  of  CCD; 
and  pixel  1500,  close  to  OH  lines)  and  for  total  power  mea¬ 
surements.  Only  observations  of  the  cold  reference  black- 
body  during  a  flight  measurement  have  been  used  for  the  cal¬ 
culation.  The  given  time  corresponds  to  the  number  of  mea¬ 
surements  multiplied  by  the  time  between  the  beginnings  of 
two  consequtivc  measurements. 
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Figure  5.  (a)  Spectrum  measured  with  THOMAS  during 
the  Second  European  Stratospheric  Arctic  and  Midlatitude. 
Experiment  (SESAME)  campaign  1995  [Titz  et  al.,  1995b]. 
Eight  frequency  channels  have  been  averaged,  (b)  Spec¬ 
trum  measured  with  the  improved  THOMAS  during  the 
CRISTA/MAHRSI  campaign  1997.  Three  frequency  chan¬ 
nels  have  been  averaged.  (The  bandwidth  of  one  frequency 
channel  corresponding  to  a  CCD  pixel  is  about  0.7  MHz.) 
The  wavenumber  scale  is  given  for  the  lower  side  band 
(LSB)  of  the  lower  panel. 


“cold  measurements”  which  is  the  minimum  requirement  for 
a  calibration  procedure  that  assumes  system  stability  (com¬ 
pare  with  section  3.9). 

3.11.  System  Sensitivity 

The  performance  of  an  ideal  submillimeter  spectrome¬ 
ter  can  be  characterized  by  the  receiver  noise  temperature 
(RNT).  Together  with  the  measurement  time  and  the  spec¬ 
tral  resolution,  the  RNT  determines  the  sensitivity  of  an  ide¬ 
al  spectrometer  [Kraus,  1966].  The  value  of  the  RNT  can 
easily  be  calculated  with  the  Y  factor  method  [Hachenberg 
and  Vowinkel,  1 982]  using  measurements  of  the  hot  and  cold 
blackbody.  The  RNT  values  achieved  by  THOMAS  during 
the  flights  have  been  around  12,000K  (DSB)  at  the  IF  chain 
output. 

The  sensitivity  of  a  real  submillimeter  spectrometer  is  de¬ 
creased  by  additional  noise  contributions  like  amplifier  fluc¬ 
tuations  or  non-Gaussian  noise  originating  in  the  first  mixer 
unit.  For  the  1997  THOMAS  measurements  the  main  ad¬ 
ditional  noise  source  was  the  first  mixer  which  caused  the 
sensitivity  to  be  about  3  times  lower  than  expected  from 
the  RNT  alone,  while  major  differences  have  been  found  for 
different  mixer  units. 

3.12.  Improvement  of  Spectra 

The  two  spectra  displayed  in  Figure  5  illustrate  the  ef¬ 
fect  of  the  THOMAS  instrumental  changes  described  above. 
Especially  the  major  increase  in  the  signal-to-noise  ratio  by 
about  a  factor  of  5  and  the  spectral  resolution  by  about  a  fac- 
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tor  of  2  are  evident.  The  shift  in  frequency  channels  of  the 
DSB  spectra  is  due  to  a  different  alignment  of  the  AOS. 

3.13.  Geometry  and  Time  Coverage 
of  Atmospheric  OH  Measurements 

The  THOMAS  observations  ranged  between  44  0  and  54  0 
northern  latitude  and  between  1  °  and  28 0  eastern  longi¬ 
tude  (central  and  eastern  Europe).  The  pointing  angles  of 
the  THOMAS  uplooking  geometry  were  85°,  80°,  and  70° 
from  zenith.  The  choice  of  the  pointing  angle  is  a  tradeoff 
between  path  length  through  the  lowermost  layers  that  are 
causing  signal  attenuation  due  to  water  vapor  and  path  length 
through  the  middle  and  upper  stratosphere  and  mesosphere 
contributing  to  the  OH  signal  [Miller  etal.,  1992],  Therefore 
a  lower  pointing  angle  (e.g.,  70°)  is  more  favorable  in  case  of 
high  water  vapor  concentrations  above  the  observer  height, 
that  is,  for  (light  altitudes  below  or  inside  the  tropopause, 
and  a  greater  pointing  angle  (e.g.,  80°)  is  more  favorable  for 
flight  altitudes  in  the  stratosphere. 

There  have  been  five  measurement  flights  on  5  days.  They 
covered  local  solar  times  (LST)  between  0700  and  2000. 

4.  Retrieval  Approach  and  Error  Assessment 

In  general,  the  nonlinear  equation  (3)  describes  the  map¬ 
ping  of  atmospheric  concentration  profiles  into  the  spectrum 
domain  as  described  in  section  2.2.  The  inverse  problem  is 
the  retrieval  of  atmospheric  concentration  profiles  from  the 
measured  spectrum,  which  is  in  general  an  ill-posed  prob¬ 
lem;  that  is,  the  solution  is  sensitive  to  small  perturbations 
in  the  measurement  [ Rodgers ,  1976].  Linearization  and  dis¬ 
cretization  yields  an  ill-conditioned  least  squares  problem. 
In  order  to  get  a  physically  meaningful  solution,  additional 
constraints  have  to  be  introduced  in  the  least  squares  prob¬ 
lem. 

There  are  several  unknowns  to  be  retrieved  from  the  mea¬ 
sured  spectrum:  First  of  all,  the  OH  concentration  profile. 
The  number  density  has  been  shown  to  be  preferable  for 
retrieval  instead  of  (he  volume  mixing  ratio  because  of  the 
smaller  altitude  dependence  causing  a  better  conditioning  of 
the  least  squares  problem  [Schimpf  and  Schreier,  1997].  In 
addition,  the  H2O  volume  mixing  ratio  has  to  be  determined 
because  of  its  influence  on  the  damping  of  the  OH  signal.  In 
contrast  to  OH,  the  volume  mixing  ratio  of  H2O  is  prefered 
for  retrieval  instead  of  the  concentration  because  of  its  lower 
altitude  dependency. 

Another  parameter  is  the  frequency  of  the  first  LO  spec¬ 
ifying  the  upper  and  lower  sideband  frequency  scale  with 
respect  to  each  other.  As  mentioned  in  section  3.4,  the  LO 
frequency  is  subject  to  drifts,  so  that  a  mean  value  is  deter¬ 
mined  by  the  retrieval  algorithm.  Last  but  not  least,  a  pos¬ 
sible  wavenumber  shift  of  the  whole  spectrum  due  to  AOS 
drift/misalignment  has  to  be  determined. 

For  the  inversions  presented  here,  the  atmospheric  pres¬ 
sure  and  temperature  input  data  were  derived  from  prelimi¬ 
nary  CRISTA  data  (version  P05)  (J.  Oberheide  and  M.  Riese, 
personal  communication,  1998)  measured  in  the  same  time 
period. 


4.1.  Retrieval  Algorithm 

Phillips-Tikhonov  regularization  using  the  L  curve  crite¬ 
rion  for  the  automatic  determination  of  the  regularization  pa¬ 
rameter  [ Hansen ,  1992]  is  a  proper  approach  to  solve  the  in¬ 
verse  problem  in  a  numerically  robust  and  efficient  way,  as 
discussed  by  Schimpf  and  Schreier  [  1 997]  in  an  atmospheric 
remote  sensing  context.  The  basic  minimization  condition  of 
the  linear  Phillips-Tikhonov  regularization  scheme  is  given 
by 

min(||Kf-g||2  +  A2||Lf||2),  (5) 

f 

where  g  is  the  measurement,  K  is  the  weighting  function 
or  Jacobi  matrix,  A  is  the  regularization  parameter,  L  is  the 
regularization  operator,  and  f  is  the  unknown  solution  of  the 
inverse  problem.  ||Kf  —  g||2  is  called  the  residual  norm,  and 
j|Lf  ||2  is  the  constraint  norm. 

Solving  the  minimization  condition  (5)  yields  a  solution 
depending  on  A 

f(A)  =K+g  =  (KtK  +  A2LTL)-1KTg.  (6) 

In  order  to  avoid  numerical  instabilities  in  (6)  due  to  explicit 
matrix  products,  minimization  condition  (5)  can  be  solved 
in  a  numerical  robust  way  using  generalized  singular  value 
decomposition  (GSVD)  [Van  Loan,  1976;  Anderson  el  al., 
1995], 

A  variety  of  methods  has  been  suggested  to  determine  the 
optimum  regularization  parameter  A,  for  example,  a  priori 
(based  upon  experience  from  previous  retrievals)  or  a  poste¬ 
riori  choice  (trial  and  error,  e.g.,  visual  inspection  of  the  so¬ 
lution),  discrepancy  principle,  generalized  cross  validation, 
or  the  L  curve  criterion.  The  L  curve  is  a  double  logarithmic 
plot  of  the  residual  versus  the  constraint  norm  depending  on 
A  and  results  in  an  L-shapcd  function.  The  comer  of  the  L 
curve,  defined  by  its  maximum  curvature,  yields  the  optimal 
value  of  A. 

L  is  the  martix  representation  of  an  operator  introducing  a 
smoothness  constraint  in  the  retrieval  problem.  Typical  op¬ 
erators  are  the  identity,  the  first-  or  second-order  derivative, 
or  the  Twomey  operator  [Twomey,  1 977] .  The  Twomey  oper¬ 
ator  asks  for  minimal  deviation  of  the  resulting  profile  from 
its  mean  value.  The  proper  choice  of  regularization  opera¬ 
tors  L  depends  on  the  actual  problem  and  can  be  determined 
by  retrievals  of  simulated  spectra  with  different  regulariza¬ 
tion  operators.  One  criterion  for  a  convenient  operator  is  the 
sharpness  of  the  L  curve  corner.  In  regions  with  few  altitude 
information  contents  in  the  spectrum  and  where  the  solution 
is  expected  to  be  constant,  the  Twomey  operator  is  suitable 
for  the  retrieval. 

In  general,  the  solution  f(A)  of  an  ill-posed  inverse  prob¬ 
lem  has  a  limited  altitude  resolution,  which  can  be  described 
by  an  averaging  of  the  true  atmospheric  profile  £ 

f(A)  =  Aft  =  K|Kft.  (7) 

The  matrix  A  is  called  the  averaging  kernel  and  plays  a  crit¬ 
ical  role  in  the  intercomparison  of  model  profiles  and  re¬ 
trieved  profiles. 
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In  the  discretized  inverse  problem  the  values  of  the  con¬ 
centration  profile  Foh  and  the  volume  mixing  ratio  of  water 
vapor  fHjO  are  given  depending  on  altitude.  Together  with 
the  parameters  of  the  frequency  shifts  of  the  first  LO  and 
the  whole  spectrum,  one  state  vector  of  unknowns  f  can  be 
composed.  For  I  representing  the  measured  spectrum,  1(f) 
being  the  discretized  fit  model,  and  K  =  being  the 

Jacobian  matrix,  (5)  can  be  regarded  as  one  Gau6-Newton 
type  iteration  step  of  the  nonlinear  problem  with 


=  w  ( 

f  Lj  05hifi  ^^8  hi  ft 

*OH 

fT 

rii2o 

r 

(8) 

=  (  K?n  K' 

Kqh 

KcH2o  )W 

_l(9) 

-  0 

0  cohLoii 

0  0 

0  ) 

ch2oLh2o  J 

W-1 

(10) 

=  W) 

-I  +  Kf, 

(11) 

where  the  superscript  c  indicates  the  current  state  in  the  it¬ 
eration  step.  Because  of  the  different  order  of  magnitude  of 
the  fit  parameters,  a  diagonal  weighting  matrix  W  was  in¬ 
troduced  to  get  a  better  conditioning  of  the  inverse  problem. 

The  parameters  coh  anti  ch2o  arc  scalar  values  constitu 
ing  the  weighting  between  the  constraints  described  by  Lon 
and  Lh:q.  Because  of  the  wcll-posedness  of  the  determi¬ 
nation  of  the  wavenumber  shifts,  these  parameters  do  not 
have  to  be  regularized,  and  the  corresponding  columns  of 
the  matrix  L  are  zero.  For  Loh  the  identity  matrix  was 
chosen  because  of  the  sharpness  of  the  L  curve  corner.  Be¬ 
cause  of  the  low  sensitivity  of  the  spectrum  to  water  vapor 
above  40  km  and  its  expected  constant  values  in  this  region, 
the  Twomey-operator  is  adequate  for  Lh2o-  The  parameters 
coh  and  ch2o  ate  given  a  priori  and  have  been  determined 
by  retrievals  on  simulated  spectra  with  different  coh/ch20 
ratios. 

4.2.  Retrieval  Error  Assessment 

Several  perturbations  occur  in  the  measured  spectrum 
which  are  only  known  with  limited  accuracy.  Most  of  these 
perturbations  influence  the  solution  of  the  inverse  problem; 
that  is,  they  result  in  errors  of  the  solution.  Summing  over 
all  perturbations  yields  the  total  error  of  the  retrieved  param¬ 
eters. 

Depending  on  their  nature,  the  perturbations  can  be  divid¬ 
ed  in  statistical  and  systematic  error  sources;  for  example, 
noise  in  the  spectrum  typically  is  a  statistical  error,  whereas 
a  base  line  offset  is  a  systematic  error.  Except  for  the  spectral 
noise,  all  perturbations  can  be  treated  as  systematic  errors. 

The  statistical  error  of  the  retrieved  profile  is  given  by  the 
projection  of  the  statistical  error  covariance  matrix  of  the 
measurement  vector  g  to  the  profile  space.  Under  the  as¬ 
sumption  of  statistically  independent,  Gaussian  distributed 
noise  of  the  measurements  g  with  a  statistical  error  covari¬ 
ance  matrix  proportional  to  the  identity  matrix,  the  error  co- 
variance  matrix  is  given  by 

C(f)=a2K+K+T,  (12) 

where  a2  is  the  variance  of  the  measurement  error. 


For  the  estimation  of  systematic  errors,  retrievals  on  sim¬ 
ulated  spectra  with  and  without  perturbation  can  be  per¬ 
formed,  and  the  difference  can  be  defined  to  be  the  system¬ 
atic  error.  The  total  systematic  error  can  be  defined  to  be  the 
square  root  of  the  sum  of  squares  of  all  systematic  errors. 

The  error  assessment  was  performed  for  an  observer  al¬ 
titude  of  12  km  and  an  elevation  angle  of  10°.  Pressure 
and  temperature  were  taken  form  CRISTA  data.  Profiles  of 
ozone  and  water  vapor  were  taken  from  the  AFGL  U.S. 
standard  atmosphere.  Because  earlier  results  of  the  MAHRSI 
experiment  showed  a  significantly  lower  OH  amount  in  the 
mesosphere  compared  to  the  AFGL  OH  profile,  the  meso¬ 
spheric  OH  concentration  of  the  AFGL  profile  was  reduced 
by  50  %  for  the  simulation.  Because  water  vapor,  ozone,  and 
continuum  emissions  are  dominant  in  the  observed  spectral 
region,  no  other  gases  were  considered. 

The  error  sources  in  the  THOMAS  results  can  be  grouped 
in  measurement  noise,  retrieval,  atmospheric,  spectroscop¬ 
ic,  and  instrument  errors  (Table  1).  Retrieval  errors  are  the 
uncertainty  in  the  ratio  coh/ch2o  and  the  uncertainty  in  the 
determination  of  the  optimal  regularization  parameter.  In¬ 
strument  errors  are  dominated  by  the  noise  of  the  first  mix¬ 
er,  a  possible  base  line  offset,  and  uncertainties  in  the  win¬ 
dow  transmission  due  to  photon-induced  absorption  (com¬ 
pare  with  section  3.2).  Spectroscopic  errors  are  uncertain¬ 
ties  in  the  spectroscopic  database  like  linestrength  and  line 
width  errors.  Atmospheric  errors  arc  uncertainties  in  the  at¬ 
mospheric  temperature,  pressure,  and  ozone  profiles.  The 
pointing  error  is  the  uncertainty  in  the  pointing  angle  adjust¬ 
ment.  Pointing  jitter  as  well  as  errors  in  the  field  of  view  of 
the  instrument  are  assumed  to  be  negligible. 

Errors  resulting  from  base  line  offset  are  compensated  by 
the  water  amount  and  result  in  wrong  damping  of  the  OH 
signal.  Errors  changing  the  effective  line  width  of  the  OH 
triplet,  that  is,  the  ILS  or  the  OH  pressure  broadening  pa¬ 
rameter  error,  yield  an  altitude  shift  of  the  OH  profile. 

Dominant  instrumental  error  sources  for  the  OH  column 
are  the  measurement  noise  and  errors  influencing  the  damp¬ 
ing  of  the  OH  signal  like  the  HjO  concentration  at  low  al¬ 
titudes,  a  baseline  offset,  the  sideband  ratio,  and  the  AOS 
nonlinearity.  Furthermore,  a  pointing  bias  due  to  the  point¬ 
ing  mirror  adjustment  results  in  a  baseline  offset  like  dis¬ 
tortion.  The  error  in  the  determination  of  the  regularization 
parameter  is  mainly  due  to  systematic  features  in  the  residu¬ 
als,  caused  by  nonwhite  noise  components  in  the  THOMAS 
spectra. 

The  expected  overall  error  budget  amounts  to  14%  and 
15  %  for  the  40-90 km  and  50-90  km  weighted  OH  column 
densties,  respectively,  whereas  24  %  are  expected  for  the  30- 
50  km  column  densities. 

5.  Inversion  Results 

The  inversion  of  the  measured  far-infrared  spectra  yields 
robust  OH  concentrations  between  about  30  km  and  90  km. 
The  resulting  profiles  can  be  considered  as  the  true  OH  al¬ 
titude  profiles  smoothed  by  the  averaging  kernel  (compare 
with  equation  (7)).  To  illustrate  the  smoothing,  Figure  6 
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Figure  6.  (left)  OH  model  profile  results  using  the  stan¬ 
dard  (model  A)  and  a  modified  chemistry  (model  B).  (mid¬ 
dle)  A  typical  set  of  averaging  kernel  rows  is  depicted.  The 
rows  are  shown  every  10  km  for  clarity.  The  averaging  ker 
nel  rows  represent  the  mapping  of  the  high-resolution  mod¬ 
el  profiles  to  the  low-resolution  domain  of  the  THOMAS 
retrieval  results.  Each  row  corresponds  to  a  given  altitude 
and  represents  the  contributions  of  the  true  OH  profile  to  the 
value  in  the  inverted  profile  for  the  corresponding  altitude. 
The  rows  are  assigned  to  15  altitude  layers  of  equal  spac¬ 
ing  between  15  km  and  85  km.  (right)  Both  model  profiles 
smoothed  with  the  averaging  kernel  so  that  they  can  be  com¬ 
pared  directly  to  the  corresponding  profile  retrieved  from  a 
THOMAS  measurement.  (Error  bars  illustrate  the  statistical 
error  only.)  The  shaded  area  represents  the  weighted  column 
density  between  40  km  and  90  km. 


(left  panel)  shows  two  model  profiles  (1400  LST,  August  7, 
1997),  a  typical  set  of  averaging  kernel  rows  (middle  pan¬ 
el),  and  the  model  profiles  smoothed  by  the  averaging  kernel 
(right  panel).  In  addition  to  the  smoothed  model  profiles, 
an  observed  profile  (1313  LST,  August  12,  1997)  is  also 
given.  As  already  mentioned  in  section  2.2,  no  altitude  in¬ 
formation  is  contained  in  the  measured  spectra  for  altitudes 
greater  than  about  50  km,  which  is  also  evident  considering 
the  shape  of  the  averaging  kernel  rows.  In  the  middle  and 
upper  stratosphere  the  averaging  kernel  rows  or  smoothing 
functions  peak  at  their  corresponding  altitude,  and  their  typ¬ 
ical  width  is  of  tire  order  of  25  km  which  corresponds  to  the 
altitude  resolution.  In  Figure  6  it  can  also  be  seen  that  the 
norm  of  the  averaging  kernel  rows  is  rapidly  decreasing  with 
decreasing  altitude  in  the  lower  stratosphere.  This  behavior 
indicates  that  the  information  content  in  the  spectra  about 
OH  in  the  lower  stratosphere  is  also  rapidly  decreasing  in 
this  altitude  region.  This  is  predominantly  due  to  the  fact 
that  the  lower  stratospheric  OH  contributes  only  very  little 
to  the  detected  spectral  radiance  compared  to  the  noise  in 
the  spectrum. 

6.  Comparison  With  Photochemical  Model 
Results 

In  this  study  we  use  the  Naval  Research  Laboratory 
(NRL)  photochemical  model  CHEM1D  to  calculate  the 
mesospheric  and  upper  stratospheric  OH  abundance  profiles 
for  the  conditions  of  the  THOMAS  observations.  This  is  the 
same  model  used  in  studies  of  OH  data  obtained  from  the 
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Figure  7.  Comparison  of  30-50  km  weighted  OH  colum- 
n  densities,  derived  from  THOMAS  observations  and  two 
types  of  photochemical  model  calculations. 


MAHRSI  experiment  [Summers  et  al.,  1996;  1997],  but  up¬ 
dated  with  more  recent  chemical  kinetics  and  cross  sections 
[DeMore  et  al.,  1997]. 

Because  of  the  different  altitude  resolution,  the  observed 
THOMAS  OH  profiles  and  the  model  results  cannot  be  com¬ 
pared  directly.  Weighted  column  densities  have  been  cho¬ 
sen  for  the  comparison.  They  are  defined  as  the  integral  of 
the  smoothed  OH  profile  between  two  altitudes.  The  shaded 
area  in  Figure  6  shows  the  weighted  column  density  between 
40km  and  90km.  Figures  7,  8,  and  9  show  weighted  column 
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Figure  8.  Comparison  of  40-90  km  weighted  OH  colum- 
n  densities,  derived  from  THOMAS  observations  and  two 
types  of  photochemical  model  calculations. 
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7.  Conclusions 

The  improved  THOMAS  instrument  was  used  to  measure 
the  thermal  emission  of  atmospheric  OH  aboard  the  research 
aircraft  FALCON  using  uplooking  geometry.  The  instru¬ 
mental  improvements  yielded  an  increase  in  the  signal-to- 
noise  ratio  of  about  a  factor  of  5  and  about  a  factor  of  2 
in  spectral  resolution.  Weighted  OH  column  densities  be¬ 
tween  30  km  and  50  km,  40  km  and  90  km,  and  50  km  and 
90  km  were  retrieved  using  a  Phillips-Tikhonov  regulariza¬ 
tion  scheme  and  the  L  curve  criterion  for  the  automatic  de¬ 
termination  of  the  regularization  parameter.  The  overall  er¬ 
ror  of  these  column  densities  was  determined  to  be  24  %, 
14%,  and  15  %,  respectively.  Comparison  of  the  results  to 
photochemical  model  calculations  using  standard  chemistry 
and  a  50  %  reduction  in  the  rale  coefficient  for  O  +  HCH  — > 
OH  +  O2  shows  that  neither  model  is  capable  of  reproduc¬ 
ing  upper  stratospheric  and  mesospheric  OH  concentrations 
simultaneously. 


Figure  9.  Comparison  of  50-90 km  weighted  OH  colurn- 
n  densities,  derived  from  THOMAS  observations  and  two 
types  of  photochemical  model  calculations. 


densities  of  THOMAS  and  model  calculations  for  different 
altitude  intervals.  The  results  of  all  five  measurement  flights 
are  included  in  the  figures  together  with  the  statistical  error 
of  the  inversion. 

For  the  interpretation  of  the  comparison  of  measurement 
and  model,  one  has  to  bear  in  mind  the  possible  system¬ 
atic  errors  and  the  significant  altitudinal  smoothing  of  the 
retrieved  profile.  For  the  given  weighted  column  densities 
the  sensitivity  to  OH  outside  the  integrated  altitude  interval 
is  about  40%  for  the  30-50  km  column,  15%  for  the  40- 
90  km  column,  and  20%  for  the  50-90 km  weighted  column 
amount. 

The  comparison  of  weighted  OH  column  densities  derived 
from  THOMAS  observations  and  model  calculations  in  Fig¬ 
ures  8  and  9  shows  that  the  absolute  values  of  the  measured 
upper  stratospheric  and  mesospheric  weighted  OH  columns 
generally  fall  between  the  standard  (model  A)  and  modified 
(model  B)  chemistry  results;  that  is,  they  are  about  15  %  low¬ 
er  than  the  predicted  values  using  the  standard  chemistry. 
Moreover,  all  THOMAS  observations  reflect  the  shape  of 
the  modeled  diurnal  variation  well.  Within  the  experimental 
uncertainty,  no  additional  influences  on  OH  other  than  the 
LST,  for  example,  local  or  temporal  variations  in  the  H2O 
abundance,  can  be  identified.  The  upper  stratospheric  and 
mesospheric  weighted  columns  also  seem  to  corroborate  the 
1994  MAHRSI  results  that  also  found  an  overprediction  of 
OH  by  standard  model  calculations  in  the  mesosphere  [Sum¬ 
mers  etal.,  1997J. 

The  larger  scatter  of  the  stratospheric  weighted  columns 
(compare  with  Figure  7)  together  with  the  larger  presumed 
systematic  error,  makes  their  clear  interpretation  difficult. 
However,  the  diurnal  variation  as  well  as  the  absolute  val¬ 
ue  of  the  stratospheric  weighted  OH  columns  is  well  within 
the  range  of  theory. 
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